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Abstract 

We examine a class of exact solutions for the eigenvalues and eigenfunctions of 
a doubly anharmonic oscillator defined by the potential V(x) = w 1 2 * /2x 2 + Ax 4 * /4 + 
yx 6 /6, y > 0. These solutions hold provided certain constraints on the coupling 
parameters to 2 , A and y are satisfied. 
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1. Introduction 

The anharmonic oscillator we consider in this note has the potential V(x) = f uj 2 x 2 + 
-jA-'c 4 + i where w 2 , A and rj > 0 are real parameters. With E denoting the energy 
eigenvalue, the associated Schrodinger equation is 

^ + (2R — w 2 x 2 — f Ax 4 — \gx Q ) ip(x) = 0 (—00 < x < 00 ) (1.1) 

subject to the boundary conditions 


lim ip(x) = 0. 

x —>-ioo 

In [1], Flessas obtained exact solutions for the energy eigenvalues and eigenfunctions 
(for the ground state and first excited state) when certain constraints on the coupling 
parameters are satisfied. Singh et al. [4] subsequently showed how an infinite set 
of such solutions could be constructed, although they presented details only of the 
eigenfunctions with 2 and 3 nodes. Again, these solutions are obtained when a sequence 
of constraints on the coupling parameters is satisfied; a finite number of such solutions 
is engendered by the fulfillment of each set of constraints. 

In this note we consider the case of eigensolutions with up to 2N and 2 N + 1 nodes, 
where iV < 3, in more detail. 
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2. Solution scheme 


We define the quantities 





c = oj 2 + (3?y) 2 


3 £ 

16 // 


With the substitution 

■0(x) = exp [— \ax 2 — \bx 4 ] y(x ) 
in (1.1), we obtain the differential equation 

y"(x) — 2(ax + bx 3 )y'(x ) + (2E — a — cx 2 )y(x) = 0. 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


We look for polynomial solutions of (2.3) of the form 

N 

y{x) = £ A n x 2n+e (A 0 = 1; N = 0,1,2,...), (2.4) 

n=0 

where e = 0 (even solutions) or e = 1 (odd solutions) and the A n are coefficients to be 
determined. Substitution of (2.4) into the differential equation (2.3) then yields 

N r 

£ A n l (2 n + e)(2 n — 1 + e)x 2n 2 + (2 E — a — 2a(2n + e)x ,2n 

71=0 ^ 


-(c + 2b(2n + e))x 2n+2 | = 0. 

The requirement that the constant terms and the coefficient of x 2N+2 should vanish 
produces 

(1 + c)(2 + e)Ai + 2 E — a(l + 2e) = 0, 

c + 2b(2N + e) = 0. (2.5) 

Thus we have the conditions 

E = ^n(l + 2e) — ^(1 + e)(2 + e)Ai (2-6) 

and 

^=© 4 ( iC-“ 2 ) =4JV+3+2e - < 2 ' 7 > 

together with the N equations for the coefficients A n 

(2 n + 1 + e)(2 n + 2 + e)A n+ i + (2 E — a — 2a(2n + e))A n 

— (c + 2b(2n — 2 + e))A n _i = 0 (1 < n < N), (2.8) 

with A n = 0 for n > N. 

In what follows, we label the eigenvalues associated with the even solutions by E- 2 m , 
m = 0,1,2,..., where Eq denotes the ground state eigenvalue, and those associated 
with the odd solutions by E^ m +i- A similar indexing applies to the eigenfunctions i/)(x), 
which will possess either 2m or 2m + 1 nodes on the interval (— 00 , 00 ), respectively. 
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(i) The case N = 0. When N = 0, we obtain from (2.6) and (2.7) the values 

Eq = ^a, 7 = 3 and E\ = | a, 7 = 5. 

The associated eigenfunctions are 

ipo(x) = exp [—|ax 2 — -j6x 4 ] and i\) i(x ) = x exp [~\ax 2 — -j6x 4 ], 

which possess 0 and 1 node, respectively. These are the special solutions of (1.1) 
obtained by Flessas [1]. 

(ii) The case N = 1. When N = 1, we obtain from (2.5), (2.6) and (2.7) for e = 0 the 
values 

E = |a — Ai, 7 = 7 (c = —46) 
and from (2.8) the single equation 

(2 E - 5a)A 1 - c = 0. 

Substitution of the above values of E and c yields the quadratic 

A\ + 2 aAi + 26 = 0, 


with solutions 

Ai = —a± \Ja 2 + 26. 

The positive root is associated with an eigenfunction with no node and so is a 
ground-state value [2, p. 21], viz. 

Eq = |a — Va 2 + 26, V’o(^) = {1 + ( Va 2 + 26 — a)x 2 }exp [—| ax 2 — |6x 4 ]. 

The negative root corresponds to an eigenfunction with 2 nodes and represents a case 
of the first even excited state: 

E -2 = fa T \/a 2 + 26, ^(x) = {1 — (\/a 2 + 26 + a)x 2 } exp [— ^ax 2 — -j6x 4 ]. 

Both cases are associated with the value 7 = 7, which gives a constraint between the 
three parameters cu 2 , A and r]. 

When e = 1, a similar procedure shows that 

E = |a — -iA] , 7 = 9 

and _ 

^4i = —|a ± 5 \/a 2 + 66 . 

Again, the positive root corresponds to the lowest odd eigenfunction (with 1 node) and 
the negative root to the second odd eigenfunction (with 3 nodes). Thus, we have 

Ei = |a — Va 2 + 66, i/ji(x) = x{l + |(\/a 2 + 66 — a)x 2 } exp [—^ax 2 — |6x 4 ], 

^3 = + \/a 2 + 66 , ip 3 (x) = x’{l — |(\/a 2 + 66 + a)x 2 } exp [—^ax 2 — 56 .x 1 ]. 

These special solutions of (1.1) were given in [4]. 
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3. Exact solutions for the case N = 2 

When N = 2, we have from (2.5), (2.6) and (2.7) for even inodes (e = 0) 

E = — Ai, 7 = 11 (c = —86) 

together with, from (2.8), the two equations for the coefficients 

12A 2 + (2 E - ha) A x - c = 0, {‘IE - 9a)^ 2 - (c + 46)Ai = 0. 

Substitution of the above values of E and c then leads after some straightforward 
algebra to 


A\ + 6 aA\ + 8 (a 2 - 2b)Ax - 16a6 = 0 , A 2 = 


2bAx 


(3.1) 


A] + 4a 

The cubic equation for A\ can be written in its reduced form as 

X 3 +PX + Q = 0, Ax = X ~ 2a, 

where p = —4(a 2 + 46) and q = 16a6. The discriminant 

A = -4 p 3 - 27 q 2 = 256{(a 2 + 46) 3 - 27a 2 6 2 } > 0, 

so that (3.1) has three real roots. Moreover, inspection of the coefficients of (3.1) shows 
that there are two negative roots and one positive root. With 


P = {-4p)2=4 


the roots are then given by 


a 2 + 46\2 


Q = 


-27 q 2 

4 p 3 


= ab 


a 2 + 46 \ 2 


6 = i arcsin Q, 


(3.2) 


Xk = Psm{6 + ^irk) (k = 0,1,2). 

It is easy to establish that the positive value of Ax corresponds to the root xi- 

Then we have the following even exact solutions: (a) a ground-state eigenvalue and 
eigenfunction given by 

Eq = |a — Psin {6 + | 7 r), V’o(^) = {1 + Aix 2 + A 2 x 4 } exp [— ^ax 2 — 56 .x 1 ], (3.3) 

where 


Ax = P sin {9 + 1 7 r) — 2a > 0, 


26(P sin (0 + | 7 r) — 2a) 
Psin {9 + | 7 r) + 2a 


> 0 ; 


(6) a first even excited state (2 nodes) with eigenvalue and eigenfunction given by 

E ‘2 = |a — Psin 9, ^ 2 {x) = {1 — |7li|x 2 — |A 2 |x 4 }exp [— ^ax 2 — jbx 4 ], (3.4) 

where 

a n ■ n o A , 2b(Psin9 — 2a) 

Ax = P sm 9 - 2a < 0, A 2 = — - . —- < 0; 

P sm 9 + 2a 
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and (c) a second even excited state (4 nodes) with eigenvalue and eigenfunction given 
by 


E 4 = |o — Psin {9 + |7r), ipi(x) = {1 — l-Ailar + A 2 X 4 }exp [—^aar — |6x 4 ], (3.5) 
where 1 


A\ = P sin (9 + 17r) — 2a < 0, 


26(Psin (0 + |7r) — 2a) 
2 Psin (0 + |7r) + 2a 


> 0 . 


A similar treatment for odd eigensolutions (e = 1) yields 
E = %a-3A U 7 = 13 (c =-106) 


and 

(3Ai) 3 + 6a(3Ai) 2 + 8(a 2 - 46)(3Ai) - 48a6 = 0, 
The reduced cubic is 


Ao = 


2bA 1 

3Ai + 4a' 


x +px + q = 0, 


3Ai = x - 2a, 


where p = —4(a 2 + 86) and q = 16a6. The three real roots are given by 


Xk = P sin (9 + ink) (k = 0,1,2), 


where now 


P = 4 


a 2 + 86 \ 2 


Q = ab 


a 2 + 86\ 2 


, 9 = 1 arcsin Q, 


Then we have the following odd exact solutions: (a) a first excited eigenvalue and 
eigenfunction (with 1 node) given by 

E\ = — P sin (9 + |7r), ip i(x) = x{l + Aix 2 + A 2 X' 1 } exp [— ^ax 2 — 46a; 4 ], (3.6) 


where 


. P . o . 2a 
^1 = j sm (9 + In) - — > 0, 


26(Psin (0 + |7r) — 2a) 
3(P sin (0 + §7r) + 2a) 


> 0 ; 


(6) a second excited state (3 nodes) with eigenvalue and eigenfunction given by 


E 3 = |a — Psin 6*, ips{x) = x{l — |Ai|x 2 — |A 2 |x 4 } exp [— |ax 2 — |6x 4 ], (3.7) 


where 


1 P ■ n 2a A 

A x = - sm 0 - y < 0, 


. 26(Psin0 —2a) 

A 2 = „ . „ / < 0; 


3(P sin 0 + 2a) 

and (c) a third excited state (5 nodes) with eigenvalue and eigenfunction given by 


P 5 = I)a — Psin ( 9 + r), ipb{x) = x{l — |Ai|x 2 + A 2 X 4 }exp [—^ax 2 — \bx A ], (3.8) 


x The quantity Psin (9 + ^n) + 2a < 0, since P| sin(#+|7r)| > (Aa/\J?>)\ sin |-7r| = 2 a when 0 < 6 < 
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where 


A\ = 


j sin (9 + 17r) 



26(Psin (9 + |7r) — 2a) 
3 (P sin {9 + |7r) + 2a) 


4. Exact solutions for the case N = 3 

When N = 3, we have from (2.5), (2.6) and (2.7) for even modes (e = 0) 

E = \a — Ai, 7 = 15 (c = —126) 

and from ( 2 . 8 ) the equations for the coefficients 

12A 2 + ( 2 E - 5a)A 1 - c = 0 
30A 3 + (2 E - 9 a)A 2 - (c + 4b) Ax = 0 
(2£-13a)A 3 - (c + 86 )A 2 = 0 . 

This produces the quartic equation for w = A\ given by 

uA + 12 aw 3 + 4(lla 2 - 156)u; 2 + 24a(2a 2 - llb)w - 366(4a 2 - 56) = 0, (4.1) 

_ 4bA\(A\ + 6 a) _ 2bA 2 

2 {A\ + 4a) (A] + 6 a) — 306 ’ 3 A\ + 6 a 

The reduced quartic is 

X 4 + PX 2 + QX + r = 0, Ai = x - 3a, (4.2) 

where p = —10(a 2 + 66 ), q = 96a6 and r = 9(a 4 + 12a 2 6 + 206 2 ). This equation possesses 
four real 2 roots. Although it is possible to express the roots in algebraic form, it was 
found that the resulting expressions were too complicated to be of practical use. In 
this case, we shall content ourselves with a numerical solution of the quartic equation 
(4.1). 

If we choose, for example, A = 0.50, 7 ] = 0.03 (a = 1.25, 6 = 0.10) then, from the 
constraint 7 = 13, we have lo 2 = 0.0625. The largest root of (4.1) together with the 
corresponding values of A 2 and A 3 are 

A x = 0.264080, A 2 = 0.021656, A 3 = 0.000558. 

We note that all the coefficients are positive and so this will result in a ground-state 
eigenfunction (no node) with the eigenvalue Eq = | — Ai = 0.360920. Values of the 
other solutions of (4.1), which we label A^ m ' > (0 < m < 3), and the associated even 
eigenvalues E 2m are presented in Table 1 and the corresponding eigenfunctions are 

1p2m(x) = {1 + A n' m) X 2n j exp [-§.T 2 - ±X 4 } (0 < TO < 3). 

n —1 
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Table 1: Values of the coefficients when A = 0.50, q = 0.03 and oj 2 = 0.0625 obtained from (4.1) and 
the corresponding even eigenvalues E^n- 


771 

^(2m) 

j (2m) 

A ( 2m ) 

E 2 m 

0 

+0.264080 

+0.021656 

+0.000558 

0.360920 

1 

-1.887128 

-0.292761 

-0.010432 

2.512128 

2 

-4.899957 

+1.859948 

+0.143071 

5.524957 

3 

-8.476994 

+8.344491 

-1.708197 

9.101994 


For odd modes with e = 1 we have 

E = §a — ‘iA\ , 7 = 17 (c = —146) 

and the equations for the coefficients 

20 A -2 + (2 E - 7a)A x - (c + 26) = 0 
42A 3 + ( 2 E - lla)A 2 - (c + 66 )Ai = 0 
(2 E - 15a) A 3 - (c + 106)^2 = 0. 

This produces the quartic equation for w = 2>A\ given by 

w 4 + 12 aw 3 + (44 a 2 - 1006)u; 2 + 24a(2a 2 - 21b)w - 1086(4a 2 - 76) = 0, (4.3) 

^ 46Ai (3Ai + 6 a) 2 b A 2 

2 ~ (3Ai + 4a)(3y4! + 6 a) - 426’ 3 ~~ 3Ai + 6 a ' 

This equation also has four real roots and, for A = 0.50, q = 0.03 and u 2 = 0.0625, the 
four sets of values of the coefficients An' n+1 ^ (0 < m < 3) and the corresponding odd 
eigenvalues E 2m+ \ are presented in Table 2. The associated eigenfunctions are 

V> 2 m+ i(x) = x{l + A^ 2 m+l ^x 2n } exp [-§x 2 - ±x 4 } (0 <m< 3). 

n —1 


5. Summary 

We have examined in detail the cases N = 2 and IV = 3 of exact solutions of the form 

N 

ij>(x) = A n x 2n+e exp [— \ax 2 — jbx 4 ] (e = 0,1) (5-1) 

Tl— 1 


2 This follows from the fact that p < 0, D = 64r — 16p 2 = —1024(a 4 + 12a 2 ?) + 456 2 ) < 0 and 
the discriminant A defined in [3, Eq. (1.11.17)] is A = 3 2 x 2 16 (a 12 + 36a 10 ?) + 402a 8 ?) 2 + 1848a 6 ?) 3 + 
3897a 4 6 4 + 35100a 2 ?) 5 + 405006 6 ) > 0. 
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Table 2: Values of the coefficients when A = 0.50, rj = 0.03 and uo 2 = 0.0625 obtained from (4.3) and 
the corresponding odd eigenvalues FXm+i ■ 


777 

^(2m+l) 

A (2m+l) 

=2 

4 ( 2771 + 1 ) 

Ebm+l 

0 

+0.243487 

+0.018657 

+0.000453 

1.144540 

1 

-0.611015 

-0.100752 

-0.003556 

3.708044 

2 

-1.690968 

+0.375069 

+0.030907 

6.947903 

3 

-2.941504 

+1.800358 

-0.271852 

10.699513 


of the Schrodinger equation (1.1), where a = A(3/?y)2/4 and b = (77/3)2 with 7 / > 0. 
The cases N = 0 and N = 1 have been given earlier in [1] and [4], respectively. It is 
found that the solution (5.1) can only exist if a constraint on the coupling parameters 
cj 2 , A and 77 is satisfied, viz. 


7 = 4 N + 3 + 2e 




For a given N and parity e of the eigenfunction, two parameters are free to be chosen 
with the third then fixed by the above constraint. 

For each value of N considered it is found that N + 1 eigenstates are produced with 
eigenvalues Eo, E?, ..., E- 2 n in the case of even modes and E\, E3, ..,, E2N+1 in the 
case of odd modes. We present below a summary of the ground-state eigenvalues and 
eigenfunctions of type (5.1) with e = 0 and normalised such that Aq = 1: 


N = 0 : Eq = io, 7 = 3, 

N = 1 : Eq = |a — V a 2 + 2b, 7 = 7; A\ = \/a 2 + 2b — a, 

N = 2 : Eq = |a—Psin(6 ) +|7r), 7 = 11; A\ = P sin(0+|7r)—2a, A 2 = ——j-^ , 

N = 3: -F 0 = |a — x*, 7=15, 

where P and 8 are defined in (3.2) and x* denotes the largest root of the quartic (4.2). 
In the case N = 4 only numerical solutions for the coefficients A n (1 < n < 3) are 
obtained for a specific choice of parameters. 


References 

[1] G. P. Flessas, Exact solutions for a doubly anharmonic oscillator, Phys. Lett. 72 A (1979) 289-290. 

[2] G. P. Flessas, On the three-dimensional anharmonic oscillator, Phys. Lett. 78A (1980) 19-21. 

[3] F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark (eds.), NIST Handbook of Mathe¬ 
matical Functions, Cambridge University Press, Cambridge, 2010. 

[4] V. Singh, A. Rampal, S. N. Biswas and K. Datta, A class of exact solutions for doubly anharmonic 
oscillators, Lett. Math. Phys. 4 (1980) 131 134. 



